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Abstract. We use the Type la Supernova gold sample to constrain the parameters of dark 
energy models namelly the Cardassian, Dvali-Turner (DT) and generalized Chaplygin gas 
(GCG) models. In our best fit analysis for these dark energy proposals we consider flat and the 
non-flat priors. For all models, we find that relaxing the flatness condition implies that data 
favors a positive curvature; moreover, the GCG model is nearly flat, as required by Cosmic 
Microwave Background (CMB) observations. 


1. Introduction 

Various proposals have been put forward to explain recent observations indicating that the 
universe is accelerating. A possible explanation is that the universe is filled with dark energy in 
the form of an exotic component, the generalized Chaplygin gas, with negative equation of state 
II1I21IHI- The striking feature of this model is that it allows for an unification of dark energy 
and dark matter [Hj. 

Another possible explanation for the accelerated expansion of the Universe could be the 
infrared modification of gravity one should expect from extra dimensional physics, which would 
lead to a modification of the effective Friedmann equation at late times. A concrete model has 
been suggested by Dvali, Gabadadze and Porrati |Hj and later generalized by Dvali and Turner 
1^. Another possibility is the modification of the Friedmann equation by the introduction of an 
additional nonlinear term proportional to p”, the so-called Cardassian model [£]. 

Currently type la supernovae (SNe la) observations provide the most direct way to probe the 
dark energy component at low to medium redshifts. Recently, supernovae data has been analysed 
by various groups and it was shown that it yields relevant constraints on some cosmological 
parameters. In particular, it is possible to conclude that, when one considers the full supernova 
data set, the decelerating model is ruled out with a significant confidence level [7j. It is also 
shown that one can measure the current value of the dark energy equation of state with higher 
accuracy and the data prefers the phantom kind of equation of state, wx < — 1. Furthermore, 
the most significant result of that analysis is that, without a flat prior, supernovae data does not 
favour a flat ACDM model at least up to 68% confidence level, which is consistent with other 
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Table 1. Best fit parameters for the Cardassian, DT and GCG models, considering flat and 
non-flat priors. 


Cardassian model 


n 



Flat Prior 

0.49 

-1.4 

— 

173.7 

Non-Flat Prior 

0.21 

-3.1 

0.47 

173.2 

DT model 




x"" 

Flat Prior 

0.51 

-19.2 

— 

174.7 

Non-Flat Prior 

0.24 

-60.0 

0.43 

174.0 

GCG model 


a 

i^k 

X^ 

Flat Prior 

0.93 

2.8 

— 

174.2 

Non-Flat Prior 

0.97 

4.0 

0.02 

174.5 


cosmological observations. In what concerns the equation of state of the dark energy component, 
it has been shown in Ref. [H] , using the same set of supernovae data, that the best fit equation 
of state of dark energy evolves rapidly from wx — 0 in the past to wx ~ — 1 at present, which 
suggests that a time varying dark energy fits the data better than the ACDM model. 

In this paper, we analyze the Cardassian, the DT and the GCG models using so called gold 
sample SNe la compilation of data by Riess et al. [2]; we consider both flat and non-flat priors. 
For more details see Ref. |1 Dj . 


2. Observational constraints from supernovae data 

The observations of supernovae measure essentially the apparent magnitude m, which is related 
to the (dimensionless) luminosity distance Dl by 


where 


Also, 


m{z) = M+ 5 logio Dl{z) , 

(1) 
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DUz) - H,{1 + z) . 

(2) 

A4-M + 51ogio(^^pJ+25, 

(3) 


where M is the absolute magnitude which is believed to be constant for all supernovae of type 
la. For our analysis, we consider the set of supernovae data recently compiled by Riess et al. 
in] known as the gold sample. The data points in this sample are given in terms of the distance 
modulus Hohs{z) = rn{z) — Mohsi^) and the is calculated from 
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where M.' = M — M^hs is a free parameter and Ca) is the theoretical prediction for 

the dimensionless luminosity distance of a supernova at a particular distance, for a given model 
with parameters c^. The errors cr^obs('^) take into account the effects of peculiar motions. 
Minimization of the y^, with respect to A4', 0^, and the respective model parameter(s) 
leads to the results summarized in Tabled The best fit value for A4' is 43.3 for all models. 














3. Cardassian model 

We first consider the Cardassian model |H1; in this model, the universe is composed only of 
radiation and matter and the increasing expansion rate is given by 
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where Mpi = 1.22 x lO^^GeV is the 4-dimensional Planck mass, b and n are constants, and we 
have added a curvature term to the original Cardassian model. The new term dominates only 
recently, hence in order to get the recent acceleration in the expansion rate, n < 2/3 is required. 



Figure 1. Confidence contours for the 
flat Cardassian model. The solid and 
dashed lines represent the 68% and 95% 
confidence regions, respectively. 


At present, the universe is matter dominated and Eq. © can be rewritten as 
/ H \ ^ 

~ + (1 “ ; ( 6 ) 

where Ho is the present value of the Hubble constant and is the present curvature 

parameter. Notice that the case n = 0 corresponds to the ACDM model. 



n n 


Figure 2. As for Figure ^ for the non-flat Cardassian model. 











For the flat case, we have only two parameters and their best fit values are = 

{0.49,—1.4}. In Fig. Qwe show the 68% and 95% confidence contour plots, where it is clear 
that the ACDM model is excluded at a 95% confidence level. 

If we relax the flat prior, we find that the best fit values are fife = 0.47, Qm = 0.21 and 
n = —3.1. Moreover (see Fig. n the case n = 0 is again excluded at 95% confidence level. 



Figure 3. Confidence contours for 
the flat DT model. The solid and 
dashed lines represent the 68% and 95% 
confidence regions, respectively. 


4. Dvali-Turner model 

The second model we consider is the one proposed by Dvali & Turner [^, where the Friedmann 
equation is modified by the addition of a new term which may arise in theories with extra 
dimensions [H] 
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where Tc is a crossover scale which sets the scale beyond which the laws of the 4-dimensional 
gravity breakdown and become 5-dimensional and we have introduced a curvature term. 




Figure 4. Confidence contours for the non-flat DT model. 


If we assume that the universe is matter dominated the Friedmann expansion law can be 




















written as 
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Notice that /3 is the only parameter of the model: for /3 


0 the new term behaves like a 





Figure 5. Confidence contours for the 
flat GCG model. 


cosmological constant and the case (3 = 1 corresponds to the Dvali-Gabadadze-Porrati (DGP) 
model [S] . The requirement that the new term does not interfere with the formation of large-scale 
structure leads to the bound (5 <1. 




Figure 6. Gonfidence contours for the non-flat GCG model. 

From Fig. |21 it is clear that j3 is very weakly constrained and can become arbitrarily large 
and negative. Moreover, both the AGDM model (/3 = 0) and the DGP model {(5 = 1), are 
strongly disfavoured. Fixing (3 = 1 (DGP model), we get Dm = 0.17 . 

If we relax the flat prior, we find that the best fit results are for {Dm,Dfc,/3} = 
{0.24,0.43,-60,0}. As we can see from the contour plots shown in Fig. @1 (3 can become 
arbitrarily large and negative, if we allow Dm to be large. Moreover, both AGDM and DGP 
models are disfavoured at 95% C.L. 







5. Generalized Chaplygin gas model 

Finally, we consider the generalized Chaplygin gas model, which is characterized by the equation 
of state 

A 

Pch — 1 (9) 

Pch 

where A and a are positive constants. For a = 1, the equation of state is reduced to the 
Chaplygin gas scenario pp. Integrating the energy conservation equation with the equation of 
state ©, one gets 0 

(1 
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a3(i+«) 

where pcho is the present energy density of GCG and ^4^ = Hence, pch, interpolates 

between a dust dominated phase, pch oc a~^, in the past and a de-Sitter phase, pch = —Pch-, at 
late times. This property makes the GCG model an interesting candidate for the unification of 
dark matter and dark energy. Notice that a = 0 corresponds to the AGDM model. 

The Friedmann equation for a non-flat unified GGG model is given by 
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This model has been thoroughly scrutinized from the observational point of view; for a GMB 
power spectrum compatibility analysis see e.g. Ref. m and for a previous supernovae data 
analysis see Ref. m and references therein. The 68% and 95% confidence level contours for the 
flat case are shown in Fig. |SI Notice that the AGDM model is consistent at 95% G.L although 
it is ruled out at 68% G.L. 

When we relax the condition of flat prior, the best fit model becomes very close to flat 
(Hfc = 0.02) which is a a new result. The confidence contours for this case are shown in Fig. 6. 


6. Conclusions 

We performed a likelihood analysis of the latest type la supernovae data for three models: the 
modified gravity Gardassian and Dvali-Turner models and the generalized Ghaplygin gas model 
of unification of dark energy and dark matter. We find that SNe la most recent data allows, in 
all cases, for non-trivial constraints on model parameters as summarized in Table ^ We find 
that, for all models, relaxing the flatness condition implies that data favors a positive curvature 
and the GCG model is nearly flat in this case. The fact that the gold sample of supernovae data 
prefers a flat GCG model, which is consistent with CMB observations, leads us to conclude that 
the GCG is a better choice among the three alternative models that we have considered. 
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